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Abstract 

We study the non-linear evolution of the curvature perturbations during matter dom- 
inated era. We show that regardless of the origin of the primordial perturbation, the 
Bardeen potential and curvature receive sizable contributions from the classical non-linear 
evolution effects, and quantify them exactly. On the super-horizon scales we have squeezed 
peak of the bispectrum with magnitude, in terms of the local non-linear parameters of 
Bardeen curvature, 1/6 < /nl ^ 19/15, and of Bardeen potential, —1/4 < /nl ^ 7/5, 
depending on the configuration of momenta. On the sub-horizon scales the bispectrum 
show equilateral shape, and can serve as a potential probe of general relativity. 
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The ongoing observation programs have been bringing the era of precision cosmology. The 
combination of the Wilkinson Microwave Anisotropy Probe on the cosmic microwave back- 
ground (CMB) and the Sloan Digital Sky Survey on large scale structure (LSS) results show 
that the primordial perturbation is characterized by a nearly flat power spectrum with almost 
perfect Gaussianity [Ij. This is consistent with the predictions of the paradigm of inflation [2], 
which provide the most successful mechanism for the causal generation of the primordial per- 
turbation on cosmologically relevant scales. The upcoming surveys on the CMB and LSS, such 
as Planck, BigBOSS, Large Synoptic Survey Telescope and Euclid, will provide data even more 
precise enough to detect small deviations, if there are any, from what the simplest single field 
inflation model predicts. They are especially important as we can gain more information crucial 
for the primordial perturbation and the generation mechanism, thereby constraining the early 
universe and high energy physics [3]. 

Among the possible deviations, non-Gaussianity has been receiving great interest for decades [1] . 
Usually, the simplest non-Gaussianity is parametrized by a dimensionless non-linear parameter 
/nl of the Bardeen curvature $, the curvature perturbation in the zero-shear gauge (often 
termed the Newtonian gauge), as [H E] 

$ = $L + /NL$i, (1) 

where the subscript L denotes the dominant linear, Gaussian component. Since this expan- 
sion is local in the configuration space, this is called the "local" type non-Gaussianity. This 
parametrization is also widely used in the studies on the inflationary non-Gaussianity of the 
comoving curvature perturbation 7?. as |6] 

7^ = 7^i + . (2) 

5 

This is based on the linear relation between $ and TZ in the large scale limit during matter 
domination 

^ = -n. (3) 

5 

This requires, however, a great caution. Primordial non-Gaussianity is one realization of 
non-linearity of the primordial perturbation. Hence, ([2]) being based on the linear relation (j3]) is 
not guaranteed to be valid. To estimate /nl and in turn the degree of non-Gaussianity properly, 
we have to reconsider ([2]) to take into account the full second order evolution effects during 
matter domination. This is important to constrain the generation mechanism of the primordial 
perturbation. For example, it is believed that the detection of local type non-Gaussianity 
would rule out any single field inflation models (TJ |8] . But this strong conclusion is based on 
the doubtful parametrization ([2]). This possibility has been anticipated and studied in a few 
literatures [U El [TJ [HI [101 HH Il2] , but they are incomplete. As we will show, in the literature, the 
non-linear evolution effects in /nl were derived often based on the linear Sachs- Wolfe relation. 
In order to estimate /nl properly, the CMB temperature fluctuations should include secondary 
effects as well pT| [T2] . Also, in the literature, the non-linear effects dominant on smaller scales 
have been neglected. With the upcoming observations, LSS provides another powerful probe to 
study non-linearities in cosmological perturbations. Thus, to constrain non-Gaussianity using 
LSS the non-linear effects on small scales are important. 
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In this article, we present the proper second order relation between $ and TZ generated by 
non-linear evolution during matter domination valid on all scales. With the correct prescription, 
/nl in $ induces a number of additional second order terms to ([2]). Further, starting from the 
primordial curvature perturbation generated during inflation, we show that the evolution effects 
give rise to sizable contributions to $. Also we make a clear distinction between the Bardeen 
curvature $ and the Bardeen potential \E', which is important beyond linear perturbation theory. 

Our starting point is the metric with scalar-type perturbations [13] 



-(1 + 2a)dt' - 2a/3^idtdx' + a' [(1 + 2^) 6ij + 2j^ij] dx'dx^ 



(4) 



where we have assumed flat background geometry and set c = 1. We take 7 = as the 
spatial gauge (threading) condition. The comoving curvature perturbation TZ and the Bardeen 
curvature $, both gauge invariant, correspond to ip — ip'^ in the comoving gauge condition v = 
with V being the spatial component of the fluid four-vector Ui = —av^i, and (p in the zero-shear 
gauge condition /3 = 0, respectively. The reason why we define 71 in such a way is because 
we want to separate the primordial non-linear component which is usually computed with the 
metric Qij = a?e'^6ij. We further introduce Bardeen potential ^ , which is a in the zero-shear 
gauge. To linear order we have \E' = — $ but as we go to non-linear order it is important to 
distinguish between the Bardeen curvature $ and the Bardeen potential ^ . Note that to linear 
order and in the sub-horizon limit to second order \E' coincides with the Newtonian gravitational 
potential. We consider a pressureless fluid. Up to second order, assuming A = 0, the relatively 
growing mode exact solutions are [121 HH fT5] 



7^ =C - - 



3^ 3 ^0 6^1 

$ =-7^^ — Tc H — A"^ 

5 25 25 



-7^A7^ + 3A-^ (7^7^'*^■) .. 
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^ = - $ + —n' + -A-i [-7^A7^ + 3A-^ inn^'^) 

25 5 L 
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(5) 

^7^'^7^,i - A-i (7^•^A7^) . 

(6)' 
(7) 



where C = C{x) is a constant coefficient of the relatively growing solution of ip in the comoving 
gauge in the large scale limit. Compared with (E]), (E]) is the correct relation between <l> and 71 
valid to second order perturbation. It shows clearly that in the matter dominated era we have 
non-trivial contributions at second order. 

In the large scale limit 7Z is known to be conserved throughout the evolution even to second 
order [151 [IS]- ^ is constant to second order for general time- varying equation of state or 
field potential. In the conventional scenario, 7Z on the super-horizon scales is generated from 
the quantum fluctuations during inflation but it may be coming from other mechanism. For 
our current purpose, it is sufficient to note that the large scale constant of dSj) represents 
the primordial component of the comoving curvature perturbation. From below we will write 
(^ — (^2 = 7? . 

On the other hand, observations are often characterized by $ or \1' which, to linear order or 
to second order in the small scale limit, satisfy \E' = — $ with being the same as the perturbed 
Newtonian gravitational potential [17]. In the context of the CMB temperature fluctuations. 
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the observation is sensitive to the non-Gaussianity of In the matter dominated era and in 
the large scale limit, with A = and flat geometry, we have |1] 

6T 

where the subscripts O and E indicate the observed and emitted epochs, respectively. 

To translate /nl given in terms of $ in ([T]) into 7^, we use ^ and ([6]) to find, in the large 
scale limit. 
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Thus, the conventional parametrization in ([2]) based on the linear relation misses substantial 
contributions. These terms describe the evolution effects after inflation up to matter dominated 
era. 

With the proper non-linear relation between TZ and $, we now proceed to evaluate /nl 



correctly. In the Fourier space, we introduce non-local and f'^l^ as 



(27r)= 
d^qid^q2 3 ^(7^) 



(27r)3 5 



&\Qi, q2)nL{q,)nM5^^\k~ q,,) , 



(10) 
(11) 



where qf^a = 9i + 92- Here, f^-^ denotes the primordial non-Gaussianity, generated from for 
example inflation in the conventional scenario. From ([5]) and (Q, we can derive 
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„ , ^ 1 kf + kj , 3(fci-fci2)' + (fc2-fci2)' 



k 



12 



"12 



"12 



1 



21^12 



k'^ P 
2ki- h2 + StJ-A^ ■ ki2 + ■ k 

k^2 f"12 



12 



(12) 
(13) 
(14) 



Apparently, the two non-local parameters and fj^ are not the same. They are related in a 
non-trivial manner, with receiving non-linear evolution effects during matter domination. 
These effects are composed of two parts: one being independent of the horizon scale k^ = aH, 
and the other being dominant on sub-horizon scales. Both are shape dependent. If we introduce 
/jlfL^ similarly defined as in (11 01) with $ replaced by \E', by using the relations between TZ and 
presented in ([5]) and ([7]), we can find 



fS^{k„ fc2) = -fi:i>{h, k^) + g{k,, h) 



tin) 



ki2 
aH 



h{ki, hi) 



(15) 



In the large scale limit this relation is presented in Refs. [H [9]. We will address this issue 
later. Both and exhibit similar structure as we can see from (fT2l) and (JTSil . and in the 
following we will closely analyze /n*^. 
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Figure 1: Contour plot of the contribution to the shape function {kik2k^Y Bi^[ki, A^, k^) of the 
horizon independent terms (upper-left), and the horizon dependent terms (upper- right). In 
the lower panel we show the corresponding shape function of We can see that in the left 
panel the shape functions are peaked at the squeezed limit. This means some of the pure non- 
linear evolution effects during matter dominated era take place over all scales, even outside the 
horizon. Meanwhile, in the right panel, the shape functions show their maximum amplitudes 
at the equilateral limit. Since these terms are heavily suppressed on super-horizon scales, the 
causal non-linear effects on the sub-horizon scales are maximized when the three momenta are 
of more or less the same size. As we can see from (fT2|) and (fTSll . on large scales the shape 
functions of \l/ and $ are similar, and on small scales the magnitude is precisely the same but 
the sign is opposite. This reflects the second order relation $ = — \1' on small scales. We have 
set ki = lOkn for the right panel. 



We can compute fj^^ by adopting for example the cubic order action [HI El [IS] or the SN 
formalism [19\ to calculate the primordial bispectrum. It is known in the literature that for 
all inflation models where only one degree of freedom is important during seed generation, in 
the squeezed limit we have the consistency relation f^^ = 5(1 — ns)/12 with Ug being the 
spectral index of the scalar power spectrum [6l [7]. Thus, it is often claimed that the detection 
of I/nl| ^ C'(l) in the squeezed configuration will rule out all classes of single field inflationary 
models. But, as stressed before, this is based on ignoring the classical contribution we have 
studied: in f|T2|) we have shown pure classical contributions in matter dominated era which 
have far larger magnitude, so that tiny primordial f^^^ is completely subdominant compared 
with /i^L^ 

We can calculate the bispectrum of $ as 

B^ih, h2, h) = -h2)P^ik,)P^{k2) + (2 perm) . (16) 
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Given $ in (HO]) and /jlj^ in 1^, we can estimate exactly the shape and the magnitude of the 
bispectrum. In Figure [H we present the dimensionless shape function (/ci/c2A;3)^S$(fci, ^2, ^3), 
normahzed by the amphtude of the power spectrum Aq> = P,^{k) / k^"'^ , and that of apart 
from possibly negligible f^- We present two contributions separately: (left) the horizon 
independent terms, and (right) the horizon dependent terms. We can see that the former is 
peaked in the squeezed limit, while the latter in the equilateral limit. Since is explicitly 
momentum dependent, in fact the perfect local non-Gaussinity ansatz, i.e. /nl is a constant, 
does not work. This suggests that /nl is, despite of its popularity, not a good parameter 
to describe the bispectrum. Nevertheless, we can estimate /nl by comparing ( fT6|) with the 
bispectrum we can find from the perfect local ansatz ([T]). Then we find the local /nl in the 
large scale limit as, with being ignored, /nl = /nl (^i, ^2, ^3) which depends on the shape 
of the triangle. It gives /nl = 1 and 1/6 in the squeezed and equilateral limits, respectively. 
But in the folded limit /nl is also dependent on n<j, and gives 19/15 for Ug = 1. Therefore, we 
have 

1 19 

g</NL<^. (17) 

In terms of /nl\ in the same way we have /nl = 1, —1/4 and 7/5 in the squeezed, equilateral 
and folded limits, respectively. Therefore, in terms of \E' we have 



-7</l 



NL 



< 
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Thus, the detection of local /nl with the amplitude of 0{1) does not necessarily exclude single 
field inflation models. From another perspective, (fT2|) and (1151) predict that, unless primordial 
is bigger than 0{1), we should find the classical contributions to /nl as a consequence of 
general relativity in matter dominated era. 

On smaller scales relevant for LSS observations the horizon dependent terms dominate and 
we have non-trivial behaviour of fj^ and f^- In Figure [2], we present the total shape functions 
for (left) $ and (right) \E' with different size of ki: (top) ki = 2kH, (middle) ki = lOkn and 
(bottom) ki = Ibkn, respectively. As one may expect, as we probe smaller scales, the horizon 
dependent contributions become more dominant, giving rise to a larger equilateral peak. This 
suggests that we have another probe of the general relativistic effects on smaller scales, where 
for example the scale dependent bias may serve as a powerful probe of the primordial non- 
Gaussianity |20] . 

Here, we compare our result with the previous studies. In Refs. [9l [10] , on large scales the 
non-linear evolution effects contribute to /nl as 

hL = fSt^-9 + l. (19) 

This obviously leads to different shape of the bispectrum as well as magnitude. For example, 
in the squeezed limit from flTIJl) we find /nl = —1/6. This discrepancy arises because ffT^ is 
based on the linear Sachs- Wolfe relation. Indeed, if we assume that the linear component of 
([8]) holds non-linearly as 6T/T = — $/3 and apply the exact solutions ([6]) and ([7j) to find /nl, 
we can obtain (fT9!) . 

In order to correctly incorporate the full second order evolution effects we must also take 
into account other secondary effects, such as the integrated Sachs- Wolfe (ISW) effects, weak 
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lensing, and so on. If /nl is confirmed to be C(l) so that the primordial contribution f^J is 
very small, it is very important to correctly identify the non-linear, second order effects. We 
would like to stress that apart from the second order transfer function, ISW, lensing etc, the 
correct relation between TZ and $ or between TZ and should be the starting point. 

There is another issue we should comment. As we mentioned before, at linear order and 
in the sub-horizon limit at second order, we have $ = — \E', and \E' coincides with perturbed 
Newtonian gravitational potential. Thus, in these two cases, up to sign convention, practically 
it does not matter which variable we choose to define the non- linear parameter /nl. In general, 
however, at proper second order this is no longer the case, as we have shown in ([6]) and ([7]). On 
small scales, where LSS serves as a powerful probe of non-Gaussianity, we use the Poisson-like 
relation between the matter density perturbation 6 and the Bardeen curvature $, 

6{k) = M{k)<l>{k) , (20) 

where A4{k) is a combination of matter transfer function, window function and so on. To study 
the effects of non-Gaussianity, we substitute the local ansatz ([1]) into $ in ( 120|) . In Newtonian 
context, $ should be in fact — \E', the gravitational potential, but even in the context of general 
relativity (120 p is also valid in the two cases mentioned above. However, at second order where 
we can properly consider non-linearity and thus estimate non-Gaussianity, especially in the 
large scale, this is no longer the case as we have shown in this work. 

In this article, we have reconsidered the relation between TZ, $ and The widely used 
relation (|2]) is properly extended to second order as (jH]) and ([7]). Using the correct second order 
relation between TZ, $ and ^, we have explicitly clarified the relation between the non-linear 
parameters in (fT2|) and (fT5|) : TZ is generated, in the conventional scenario, during inflation 
and the the Bardeen potential \l/ is directly related to the temperature fluctuations as (|8]) and 
the Newtonian gravitational potential. The observationally relevant /nl contains substantial 
non-linear evolution effects during matter dominated era besides the primordial contribution. 
While the one effective on super-horizon scales gives the maximum amplitude of the bispectrum 
at the squeezed limit, the one dominant on sub-horizon scales gives an equilateral peak. This 
suggests, in addition to our main findings in this article, another interesting way of probing 
general relativistic effects using non-Gaussianity on different scales. 
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Figure 2: Contour plot of the full shape functions for (left) $ and (right) \1/ with (top) ki = 2kH, 
(middle) ki = lOkn and (bottom) ki = 15kH- Depending on the scale we probe, we may have 
a mixed shape of the squeezed and equilateral limits. On smaller scales, the causal, sub- 
horizon gravitational interactions become important and as a result we find a larger peak at 
the equilateral limit. The contour scaling is the same as Figured] 
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